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ABSTRACT 



A digital simulation language for the interactive definition and 
solution of piece-wise continuous non-linear ordinary differential 
equations using the Runge-Kutta method has been designed and imple- 
mented. The combination of interactive graphics approach and a 
special differential equation description language make the analysis 
program very versatile and easy to use. For second order systems, 
a grid of phase trajectory segments over the user specified phase 
plane is used as a background for the solutions. 
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I. INTRODUCTION 



The goal of this study was to develop an interactive analysis tool 
for the study of piecewise continuous non-linear ordinary differential 
equations. It is recognized that there are available in most computer 
centers several programs for solving non-linear differential equations. 
However they are either very special purpose or else require that the 
problem and/or control description be inserted into a subroutine or 
part of the analysis program. This in turn must be compiled by the 
computer before execution time. If the user wants to change some 
feature of the system, he must change the applicable program cards 
and recompile. Thus in general only one problem type per computer 
run is possible. 

Also there exists a large number of continuous digital simulation 
languages [Clancy and Fineberg, 1965]. These are best represented 
at this time [Cardenas and Karplus, 1970] by MIMIC, DSL/90, and 
CSMP. DSL/90 and CSMP, which is based on DSL/90, have transla- 
tors which convert the prescribed input language into FORTRAN sub- 
routines. These subroutines are then compiled by the FORTRAN 
compiler and executed with the solution routine. This is a multi-step 
process which is both time consuming and not compatible with the 
interactive approach. CSMP and DSL/90 will accept the problem 
description either in the block diagram form or the differential equa- 
tion form. There is a version of CSMP (similar to the above CSMP 
in name only) which uses the block diagram approach only and is 
suitable for use on the small computers of the IBM line. MIMIC 
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incorporates its own compiler and has a very flexible algebraic 
capability. 

The program described in this thesis uses the differential equation 
approach and has a built-in compiler flexibility similar to MIMIC. In 
addition it has been designed from its conception to be fully interactive 
using graphics as the interface between the user and the program. All 
problem description is handled as data. Thus only the compiled version 
of the analysis program is needed by the user. This is accomplished 
by the use of a special differential equation description language. An 
automatic compiler (transparent to the user) built into the program 
provides the interface between the description language and the fourth 
order Runge-Kutta solver. This use of a special differential equation 
description language and the interactive graphics approach give the 
program considerable flexibility in the type of differential equations 
it can handle. 

The class of piece-wise continuous non-linear ordinary differential 
equations handled by this program is restricted by the following con- 
side rations : 

1. The equation must have only one independent variable. 

2. The equation must be normalized so that the coefficient and 
exponent of the highest order derivative are both one. 

3. At each given moment during the solution one and only one 
differential equation can be applicable. The applicable differential 
equation must be well defined for the conditions at that point. Note 
that this allows piece -wise continuous differential equations with 
discontinuities but excludes simultaneous differential equations. 
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These are the only restrictions inherent in the method of solution. 
The implementation described in this thesis has been arbitrarily 
limited to differential equations of eighth order or less. For a given 
problem there may be added restrictions in defining the differential 
equations due to the limitations of the program language used to define 
the systems and/or the ingenuity of the user. However, if the user 
can define his system using the defined programming language and the 
above three requirements are not violated then the program will solve 
the system. 

Although the language used in the system description is quite 
natural to most users, it is important to completely specify it. Thus 
Chapter III deals with a description of the elements of the language 
and the principles of its use. 

Other features of the program are: 

1. The program has the generality required to be useful for a 
wide variety of applications. In addition, the program does not re- 
quire sophisticated computer techniques on the part of the user. 

2. The outputs of the program are a smooth x vs x 1 plot (phase 
plane in the second order case) of the solution and printer plots of 
the time response of each state variable. 

3. The "window" represented by the x vs x 1 plane can be adjusted 
in size, scale, and position. These three parameters can be adjusted 
in the x and x 1 directions independently; thus the user has full flexibil- 
ity to adjust the window as desired. 

4. The solution is accomplished using the state variable approach 
with fourth order Runge-Kutta integration. The parameters of the 
solution are completely under the control of the user. 
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5. For second order time independent systems there is super- 
imposed on the phase plane a grid of phase trajectory segments. 

These slope markers represent the slope of the trajectory if it were 
to pass through that area of the phase plane. Thus without obtaining 
a solution one can visually obtain a feel for what the system will do 
for various initial conditions. This in turn aids the user in deciding 
what the essential characteristics of the system are. 

6. For the Interactive Graphics Version, the parameters of each 
feature mentioned above are available for modification as desired on 
an interactive basis. For example one may select specific values of 
various parameters and observe their effect on the resultant solution. 
The program is interactive and prompts the user in the choice of 
responses indicating his desires. 

There are two versions of this program. One version is called 
the Batch Version and runs on the IBM 360/67. The other is an 
Interactive Graphics Version and is implemented on the Electrical 
Engineering Department's XDS 9300 with the associated ADAGE 
Graphics Terminals. The interactive graphics version was primarily 
practical due to the power of the "intelligent 11 Adage Terminals. 

In Chapter II a general description of the program and its various 
parts are discussed. A major portion of the coding is the compiler 
and interpreter which are one of the primary features of this program. 

Chapter IV concerns itself with a description of the program from 
the users viewpoint. In addition to the general comments, there are 
included sections which deal separately with the Batch Version and 
the Interactive Graphics Version. 
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The use of the Interactive Graphics Version requires a knowledge 
of how to use the graphics system. This is explained by a suggested 
first session exercise. In this exercise the user is led step-by-step 
through all the techniques required to use the program. 

Chapter V ties all the above together by presenting a large number 
of actual problems that have been studied with this program. The 
input deck is explicitly presented along with the resultant output. 

A special benefit of the analysis feature of this program is its 
usefulness as a teaching or laboratory tool. This is especially true 
for the Interactive Graphics Version used to analyse a second order 
time independent system. For example the user could type in the 
description for a simple second order pendulum system. The slope 
lines give a very dramatic overall view of the effect of various types 
and amounts of damping. Since the system is interactive, the user 
can immediately pursue and verify theoretical observations. 
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- II. GENERAL DESCRIPTION OF THE PROGRAM 



A. INTRODUCTION 

The program consists of five primary parts. The first four are 
the compiler, interpreter, slope line generator (used in the second 
order time independent case only), and the solution generator. Each 
of these will be explained in more detail later in this chapter. For the 
present, a brief functional description will be given. 

The "Compiler" takes the differential equation or system of dif- 
ferential equations, as the case may be, and through a process to be 
described later comes up with a polish string [Forsythe, 1969, Bur- 
roughs, 1964] of integer codes stored as a single vector. This polish 
string is the compiled version of the input system of equations. 

The "Interpreter, " when directed, uses the polish string and the 
current value of all its required arguments and produces a value for 
the derivative of each state of the system with respect to time. The 
required arguments are indicated by the system equations and may 
include the current values of the states, time, and the values assigned 
the constant coefficients A through H. 

The resultant state variable derivatives as defined by the inter- 
preter are needed by both the slope line generator and the solution 
generator. The slope line generator calculates the grid of slope lines 
as shown in Chapter V. The solution generator, as the name implies, 
generates the time solution. 

The fifth primary part of the program is the coordinator. It is 
the coordinator that ties the above together and makes it work. It 
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provides the required scheduling based on the results of the user's 
modifications. These modifications are accomplished on an interactive 
basis when using the XDS/AGT Interactive Graphics Version and via 
data cards on the IBM 360/67 Batch Version. 



B. 



MATHEMATICAL BASIS OF THE PROGRAM 

1. The General Mathematical Concept 
th 

Assume an n order differential equation of the form 

+ f(x) = g(x,u) (2-1) 

dt n 



where f(x) is a function of the n-1 derivatives, and g(x, u) is a function 
of the inputs and the n-1 derivatives. This equation can be rewritten 
as n first order differential equations. Let us assume the following 
notation and substitutions 

x = Z(l) 

x 1 = Z(2) = ZDOT(l) 

x" = Z(3) = ZDOT(2 ) 

x'" = Z(4) = ZDOT(3 ) (2-2) 

• • 

• • 

• • 

x (n-1) = Z(n) = ZDOT(n-l). 

Rewriting equation 2-1 using the substitutions and notation specified 
in Equations 2-2 gives 

ZDOT(l) = Z(2) 

ZDOT(2) = Z(3) 
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(2-3) 



ZDOT(n-l) 

ZDOT(n) 



Z(n) 

g(x, u) - f(x). 



Each of the Z's will be called the state variables of the system. 

Note that for the first n-1 state variables the derivative of the i state 

variable of the system is equal to the i+1 state variable. The deriva- 
th 

tive of the n state variable is equal to g(x, u) - f(x). 

The solution generator needs only these derivatives of the state 

variables of the system. The compiler and interpreter use the above 

observations as the basis for the method of implementation. 

Assume that equation (2-1) has been presented to the system. 

d n x 

The first action of the compiler is to recognize the term ^n which 
is denoted on the data card as an M x M followed by n apostrophes. This 
notation is completely analagous to the use of dots to denote time de- 
rivatives. The notation will be explained in detail in Chapter III. 

Having recognized this term it now knows the order of the differential 
equation. Next the expression g(x, u) - f(x) is compiled. The compiled 
version will be used by the interpreter to evaluate the desired 
derivatives . 

The interpreter first evaluates the compiled version of g(x, u) - f(x) 
which becomes the value of the highest order state derivative. The 
other derivatives are defined using the state values as defined in 
Eq. (2-3). 

2. Mathematical Derivation of the Slope Lines 

The slope lines are used only for second order time indepen- 
dent systems. Thus the comments of this section will apply to a 
second order system only. The general system can be written as 
x 11 + f(x) = g(x) (2-4) 
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where f(x) + g(x, u) are as defined in equation 2-1. 

Letting x = Z(l) and x' = Z(2) gives 

dx/dt = ZDOT(l) = Z(2) 

dx'/dt = ZDOT(2) = g(x, u) - f(x). 

For a given time t the slope of the solution trajectory 



is 



slope = (dx/dt)/(dx‘ /dt) 
= ZDOT(2)/Z(2). 



( 2 - 6 ) 



Thus knowing a position and a corresponding slope the desired lines 
can be generated using simple trigonometry. 



C. COORDINATOR 

Figure 2-1 is a simplified block diagram of the program and out- 
lines the function of the coordinator. The coordinator provides the 
overall coordination and interface for obtaining the desired chain of 
events. The program is best summarized by reviewing this block 
diagram. 

The first action of the program is to obtain the system and plot 
description data. This includes the system equations, the plot para- 
meters such as scales and sizes, and the solution information such as 
the initial conditions, initial and final times. With this information a 
complete run through the program can be made. 

Next the compiler compiles the system description data. The 
results of this compilation, if successful, are saved in a vector (a 
1 -dimensional array or stack). If the result of the compilation or one 
of the steps which uses the compiled version encounters an illegal 
situation then an appropriate message is presented and the program 
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Fig. 2-1 OVERALL BLOCK DIAGRAM 



goes on to get new data. For the Batch Version the next new problem 
would be started. For the Interactive Version the user would make 
the appropriate corrections and proceed. 

After a successful compilation the plot is started and the grid of 
slope lines is computed and plotted. The time solution is then com- 
puted and plotted. 

At this time, the program is ready for the next set of data from 
the user. Depending on what new data is obtained several courses of 
action are taken. If the system description equations are changed 
then a recompilation is required. The program flow after recompila- 
tion is the same as before. If the system equations were not changed, 
then the coordinator checks for a change in the constant coefficients 
or the plot parameters. A change in these will cause a transfer to 
the slope generating segment. 

The Batch Version differs slightly at this point. Instead of asking 
if the "coefficients or plot parameters changed? , " the question should 
read: "Is the next solution desired on a new graph?" This allows a 

family of solutions to be plotted on the same graph if so desired. 

The above cycle of events continue until the user is done. 

D. DESCRIPTION OF THE COMPILER 

The programming details of the compiler are not of central im- 
portance to this thesis. A copy of the coding is in the overall program 
printout included as Programs A and B. Frequent use of comment 
statements including explanation of codes, key variables, and such 
are included. 

Referring to the simplified block diagram, fig. 2-2, a very brief 
description of the compiler will be given here. The compiler accepts 
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Fig. 2-2 COMPILER SIMPLIFIED BLOCK DIAGRAM 



the input a line at a time. The recognition of the elements of the 
language is accomplished by brute force. For example, let us assume 
that a n T M has been recognized. This T could be the variable T or 
the first letter of either "TAN" or "THEN". A check of the following 
characters quickly reveals which of the three it is. When a language 
element is recognized a corresponding integer code is placed in a 
stack. At the end of the recognition phase, the entire set of equations 
will have been coded into a string of integer codes stored in the stack. 

A blank card or line signals that the recognition phase has been 
completed. At this time the string is in the form known as "infix. " 

The final task of the compiler is to convert this infix string to a 
"polish" string. The polish string has the advantage that it is quite 
easy to implement an interpreter for it. This interpreter will be 
discussed in the next section. 

A part of the compilation work is to recognize and use information 
such as the order of the differential equation. The order of the dif- 
ferential equation is determined from the equations and is required by 
the interpreter. 

E. DESCRIPTION OF THE INTERPRETER 

The purpose of the Interpreter is to take the compiled code which 
was generated in the compiler and values for the required operands 
which were specified by the user and to provide the derivatives that 
are needed by the slope and solution generators. The use of these 
derivatives will be explained in the next two sections. 

The functional block diagram is presented in figure 2-3. Upon 
entry to the Interpreter all the required operands have been defined. 
These are the present values of each of the state variables, the values 
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Fig. 2-3 INTERPRETER SIMPLIFIED BLOCK DIAGRAM 
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of the constants (A through H), and the problem time T. The Inter- 
preter takes one by one each coded entry in the polish stack (result 
of the compilation phase) and carries out the required operation. The 
various types of actions are summarized by the block diagram. 

The code may be carrying information such as the order of the 
differential equation and whether the differential equation to be used 
is a function of time. These codes cause the proper variables and 
flags to be set. 

The logical branching is done in the following manner. Each 
statement in the stack is separated by a code representing a n ; n . If 
the interpreter evaluates the argument of an IF statement to be false, 
then it skips through the code to the entry following the next ,, ; n . 

Thus it has skipped the part of the statement indicating action to be 
accomplished if the logical argument is true. The SKIP instruction 
uses the same method but will skip several instructions as requested 
by the programmer. 

The arithmetic evaluation of the polish stack is accomplished by 
the use of software LIFO (last in first out) stacks. The basic prin- 
ciples involved in evaluating a polish stack are well known [Rosen 1967, 
Forsythe 1969]. For ease in implementation three different LIFO 
stacks have been used. One is used for the arithmetic operations, one 
for logical operations, and one to hold addresses. Thus they are 
respectively real, logical, and integer stacks. Associated with these 
stacks there are the required pointers to keep track of their top 
positions. 

When the interpreter notes that a differential equation has been 
evaluated it terminates the interpretation phase. As will be explained 
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in Chapters III and IV the program uses the first differential equation 
whose conditions for use are met. There it will be shown that the 
program easily handles a system whose behavior must be described 
by different differential equations in different portions of the solution 
space. The last action of the Interpreter is to define the desired 
system derivatives as described in Section B. These derivatives plus 
the order of the system are passed back to the calling slope or solution 
generator. 

F. DESCRIPTION OF THE SLOPE GENERATOR 

The Slope Generator presents a grid of slope markers over the 
area of the phase plane. The use of these slopes is demonstrated in 
Chapter V. The mathematical basis of the slopes was described in 
Section B of this chapter. 

For each grid point on the phase plane the interpreter is called. 

If the system, for that grid point, is not second order or time indepen- 
dent then it skips on to the next grid point. The effect of this skip is 
to suppress the slope markers at that point. If the system is of second 
order and not time dependent then the slope is defined and is drawn 
at the current grid point. 

The resultant derivative information from the interpreter and the 
position information are used to calculate the desired slope. The 
slope is used in turn to compute the vector that represents the slope 
at that point. This derivation is made somewhat more complicated 
by the scaling options available to the user. Thus both the slopes and 
the vector length calculation must consider the horizontal and vertical 
scale factors. 
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G. DESCRIPTION OF THE SOLUTION GENERATOR 

.The Solution Generator incorporates a very straightforward use 
of the Runge-Kutta method. The heart of the Runge-Kutta solution 
is the function subroutine RKLDEQ. This function subroutine is part 
of the IBM/360 library. For flexibility the actual cards are included 
in this program. 

After each solution point is generated the result is stored for 
future plots. In addition, the interactive graphics version presents 
the solution as it takes place. The solution rate for a typical system 
is about ten iterations per second with graphics. Thus the resultant 
presentation is approximately real time for a solution time step of 
0. 1 seconds. 
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III. GENERAL DESCRIPTION OF THE LANGUAGE 



A. INTRODUCTION TO THE LANGUAGE 

The language primarily resembles the familiar (to most) form of 
FORTRAN. Anyone familiar with one of the languages; FORTRAN, 
Algol, Basic, or PL/l will find the language used here quite easy. 

The requirements of this language are very limited and of a special 
nature. Thus, since we have a special purpose language, it is not 
complicated by the generality required of most languages. 

B. ELEMENTS OF THE LANGUAGE 
1. The Character Set 

The following characters are used in this language. 

a. The alphabetic characters as applicable. 

b. The numeric characters, 0 through 9. 

c. The special characters: 

❖ 

f 

9 

( 

) 

+ 

/ 

blanks 

d. Blanks are optional and may be used as desired for 
clarity and neatness except in the middle of a word, variable name, 
or number. The semicolon is used as a delimiter allowing more than 
one statement per line and may optionally be added at the end of a line. 
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2. Constants 



A constant is a fixed, unvarying quantity. The constants used 
in this program are of two general types. They can either be self- 
defining numeric values appearing in the source statements, or con- 
stants defined externally to the program. 

There are eight externally defined constants with names A, 

B, C, . . . H. The values of these real constants are defined by the 
user externally to the system statements. Thus the values of these 
constants may be changed as desired, between executions, without 
recompilation of the systems equations. The method of defining these 
constants is explained in Chapter IV. 

The self-defined constants are recognized by the compiler 
during compilation. With the exception of integer powers all self- 
defined constants are treated internally as real. However, when 
used in the program, the use of a decimal point is optional. That is, 
the decimal point may be either explicit or implicit. 

In the case of exponentiation, an integer, positive or negative, 
used as a power will yield the expected integer power results. On the 
other hand, non-integer powers are computed using the natural log- 
arithm/exponentiation method. 

Examples : 

Valid real numbers: 9. 05, 0. 06, 78, etc. 

Valid use of integers: 9. 7**4, V12**-3, V12 ** - 3, etc. 

Valid use of non-integer exponents: V23**0. 897, etc. 

3. Variables 

A variable is a named quantity whose value may change 
during execution of the program. Variables are specified by name. 



26 



The name must start with the letter "V" followed by one or two 
decimal digits representing the number 1 to 40 inclusive. These 
variables are always considered as real variables. 

Example: VI, V23, V39, etc. 

In addition the letter M T M is used as a variable. T is the real 
variable which represents the problem time; T is the time variable of 
the system solution. M X n is the variable of the differential equation 
and will be explained later. 

4. Expres sions 

An expression is a string consisting of constants, variables, 
and function calls interspersed with operators and parentheses. The 
formulation of an expression is governed by normal mathematical 
convention and the rules given below. There are two kinds of expres- 
sions; arithmetic and logical. The value of an arithmetic expression 
is a real number. The value of a logical expression is always a truth 
value, TRUE or FALSE. Expressions may appear in If Statements, 
Assignment Statements and in Differential Equation Statements, 
a. Arithmetic Expressions 

(1) The arithmetic operators are as follows: 

Arithmetic Operator Definition 

** Exponentiation 

£ Multiply 

/ Division 

+ Addition 

Subtraction 

(2) Rules for the Construction of Arithmetic Expressions: 
(a) All desired computations must be specified 
explicitly. For example, to multiply the variables VI and V3 one must 
write VI *V2 and not V1V2. 
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(b) No two numeric operators may appear in sequence. 

(Note: ** is considered as a single operator). Thus A* - V19 is 

illegal but may be corrected by the addition of parenthesis obtaining 

A*(-V19). 

(c) Parenthesis may be used as desired for clarification 
and to force the proper order of computations. 

(d) A mathematical function, together with its argument, 
may be used and is treated as a variable. 

(3) Order of computation in Arithmetic Expressions: 



Operation Hierarchy 

Evaluation of functions 1st 

Exponentiation 2nd 

Multiplication and division 3rd 

Addition and subtraction 4th 



(4) Examples of Arithmetic Expressions: 

(A*V3 + 5. 7)/EXP(Vl*B*T) 

-A*(l. 0 - X**2)*X ! - B*X 
A*X ! *ABS(X') + SIGN(X) 

b. Logical Expressions 

(1) Relational Operators. There are six relational operators. 
They, and the logical operators which follow, must always start and 
end with a period. The relational operators are binary operators 
which operate on arithmetic expressions and express an arithmetic 
condition which can be either true or false. The relational operators 
are as follows: 

Relational Operator Definition 

. GT. Greater than 

. GE. Greater than or equal to 
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. LT. 


Less than 


. LE. 


Less than or equal to 


. EQ. 


Equal to 


. NE. 


Not equal to 



Example s 
VI . LT. X* 

0. 94 . GE. EXP(A** - 0. 5) 

X . LT. -A .OR. X . GT. A 



(2) Logical Operators. The three logical operators are 
.NOT. , .AND. , and .OR. . Their use is defined in the 
ordinary sense of logic. 

(3 ) Rules for Construction of Logical Expressions: 

(a) Relational operators may only operate on expres- 
sions which, when evaluated, have a numeric value. 

(b) Logical operators may only operate on expressions 
which, when evaluated, have a logical value (i. e. , TRUE or FALSE). 

(c) Relational operators may not appear in sequence. 

(d) Two logical operators may appear in sequence only 
if the second one is the logical operator . NOT. 

(e) A relational and a logical operator may not appear 
side by side in an expression. 

(f) As in Arithmetic Expressions, parenthesis maybe 
added as desired for clarification and to force the proper order of 
calculations . 



(g) If the logical operator . NOT. operates on an ex- 
pression with more than one term, the expression must be enclosed 
with parenthesis. 
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(4) Order of Computations in Logical Expressions. 



Operation Hierarchy 

Evaluation of functions 1st 

Exponentiation 2nd 

Multiplication and division 3rd 

Addition and subtraction 4th 

. LT., . LE. , . EQ. , . NE. f . GT. , . GE. 5th 

. NOT. 6th 

. AND. 7th 

. OR. 8th 



If two operators have the same hierarchy, then they are used in the 
order they appear left to right. 

(5) Examples of Logical Expressions . 

(Vl**2.3 . GE. A) .AND. (X 1 . LT. 0) 

T .GT. 1.5 

T . GT. 1.5 .OR. (X .NE. 0.5) .AND. .NOT. (X . LT. X') 
(T .GT. 1.5 .OR. (X . NE. 0.5)) .AND. .NOT. (X . LT. X') 
These last two are not equivalent since the added set of parenthesis 
has altered the computation order in the last case. 

5. Mathematical Functions 

The language provides for the use of various mathematical 
functions. The present library is listed below. One calls these 
functions in an implicit manner by simply using the proper function 
name, followed by its argument in parenthesis, as an ordinary 
variable. The argument itself may be any valid mathematical ex- 
pression which in turn may contain function calls. 

a. Examples of Use. 

EXP(Vl) 

EXP(V1*V3/SQRT(SIN(X/X'))) 

(SIN(X) + COS(X**2))/3. 75 
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b. The various available functions and their preassigned 



names are expressed below. 

Name Description of Function 

SIN Sine of an angle in radians 

COS Cosine of an angle in radians 

TAN Tangent of an angle in radians 

ABS The absolute value of the argument 

EXP Exponential of the argument 

LN The natural logarithm of the argument 

LOG The base 10 logarithm of the argument 

INT Truncates the real number to an integer then reconverts 

it to a real. 

SIGN Magnitude of 1. 0 with the sign of the argument. 

6. Arithmetic Assignment Statements 

The goal in developing the above language elements was to 
make use of them as the building blocks for Arithmetic Assignment 
Statements, Control Statements, and Differential Equation Statements. 
This section will discuss the Arithmetic Assignment Statement. 

The end result of using an Arithmetic Assignment Statement 
is to compute a new value for a variable. The general form of an 
Arithmetic Assignment Statement is a = b 

where: a is a valid variable name (written without a sign). 

b is any mathematical expression following the above 
guidelines . 

= is an order to compute the value of the expression b 
on the right, and then give that value to the variable 
named on the left. 

A lengthy discussion can be made on the difference between the 
Arithmetic Assignment Statement and a mathematical statement of 
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equality. It is perhaps unfortunate that an " = " is traditionally used in 
both. cases. However if one keeps the above definition in mind, there 
should be no confusion. For someone completely new to computer 
programming methods, it is recommended that one of the many good 
fundamental FORTRAN manuals be consulted. The formation of and 
use of the Arithmetic Assignment Statement is identical with that of 
FORTRAN. 



Examples: 

1 2 

Let VI = y— + 

Zx . Z 
4x 



then the coded version is: 



VI = -1. /(Z. *X) + A**2/(4*X**2) 

2/5 

Let V5 = (A + Bx ) 



then the coded version is: 

V5 = (A + B*X**3)**(2. 0/5. 0) 
Let V9 = (SIN 2 x + COS 2 x)' 
then the coded version would be: 



V9 = (SIN(X)**2 + COS(X)**2)**0. 5. 

Let V23 = — + log I tan ^ I 

cosx ° I Z I 

then the coded version would be: 



V23 = 1. 0/COS(X) + LOG(ABS(TAN(X/2. 0))). 

Let V5 = V5 . , + | ax 2 
new old 



then a coded version would be: 



V5 = V5 + ABS(A*X*X). 
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7. Control Statements 



a. The Skip Statement 

The control statements are of two general types. One 
is a simple unconditional transfer statement called a Skip Statement. 
Its form is: 

SKIP N 

where: N is an integer that specifies how many statements should be 

skipped. 

b. The Arithmetic IF Statement 

The other more powerful control statement is the 
Arithmetic IF Statement. This is essentially the statement that one 
is familiar with in FORTRAN. Its form is: 

IF(a) THEN b; 

where: a is any valid logical expression as previously defined. 

b may be an Arithmetic Assignment Statement, a Skip 
Statement, or a Differential Equation Statement. 

It commands the computer to evaluate the logical expression a. If a 
is true then do b, otherwise skip b and go on to the next statement. 

The use of the keyword THEN is optional as is the n ; n . 

c. Examples of Valid Control Statements 

SKIP 3 

IF(X .GT. X 1 ) THEN SKIP 4 

IF(X .LE. ABS(X')) THEN Vl = A*X**2 

IF(X . GE. 0.0) THEN X !, + 0. 2*X' + X = 0. 0 

8. Differential Equation Statements 
a. Notation 

Differential Equation Statements are intended to be as 
natural as possible. The apostrophe notation, which is widely used, 
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has been adopted. Thus this program used x M as the symbol for the 
second derivative of x. 

Other derivatives of X are handled in an analogous manner. 
With the exception of this notation for the derivatives, the remainder 
of the differential equation is written using the construction methods 
used for mathematical expressions. Note that in the case of a Dif- 
ferential Equation Statement, the statement actually expresses an 
equation. Thus for this statement only, the equal sign denotes 
equality rather than an assignment. 

Example: X ,M + A^X , ^ABS(X l ) + X**2 = V1*T 

b. Writing a Differential Equation Statement 

The form of a Differential Equation Statement is: 

a — b - c 

where: a denotes the highest derivative of X written using 

the notation of the previous section, 
b and c denote any valid mathematical expressions 
using any or all of the constants and variables and 
using any or all of the lower derivatives of X as 
variables. 

It can be seen that the rules for writing b and c are extremely liberal. 
Thus the type of differential equation that can be handled is very 
general. The differential equations may be extremely non-linear, 
and have coefficients that vary with time and/or position. 

c. Examples 

The above generality can best be demonstrated by a few 
examples. Valid constructions are: 

X m + A*X M + 0. 75*X* + 0.0 
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X M + (X*X')**X + SIN(X) = 10.0 
X ,?M . + X MT + X M /X = LOG(X ,M **0. 5) 

C. USING THE LANGUAGE 

1. General Comments 

The goal of all the above elements of the language has been to 
develop tools for use in defining the differential equations of a system. 
The above language features make it possible to handle a wide variety 
of system types. This section demonstrates, by use of several ex- 
amples, how to use the elements of the language to describe various 
systems. The series of examples will start with simple applications 
and then progress into more involved descriptions. 

2. A Simple Differential Equation 

The simplest, but still very useful and powerful, form of 
system description consists of one Differential Equation Statement. 

Let us take for example Van der Pol*s 
Equation: X M - A*(l. 0 - X**2)*X' + B*X = 0. 0 

This one data card, or line on the graphics terminal, completely 
describes the differential equation with the exception of A and B. 

The variables A and B are specified on an additional data card. When 
it is desired to alter the values of A and B they are changed independ- 
ently of the system equations. Considerable execution time is saved 
by eliminating recompilation of the system equations when it is only 
desired to alter the value of a coefficient. If, however, one did not 
anticipate varying coefficients A and B, then the differential equation 
could be coded as follows, for example: 

X n - 0. 95*(1. 0 - X**2)*X* + 1.7*X =0.0 
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Assume now that it was desired that the first coefficient be 0. 90 vice 
0. 95, There is still no problem, the required recompilation is trans- 
parent to the user. The on line user will only notice a somewhat 
longer delay time before the solution starts. 

3. More than One Differential Equation Statement 

In most control system applications more than one differential 
equation is needed to describe a given system. However, only one 
differential equation is applicable at a given moment in the solution. 
For example, due to a particular type of non-linearity, the selection 
of the proper differential equation may be dependent on the sign of the 
velocity or position. In another situation, the selection of the dif- 
ferential equation may be a function of time. Numerous examples of 
various control system descriptions appear in Chapter V. 

For purposes of illustration, we shall assume the following 

system: 

If X is less than -0. 2 
then the differential equation is 

X n + 0. 2*X' + X = 0. 2 

If X is greater than 0. 2 
then the differential equation is 

X M + 0. 2*X' + X = -0. 2 

If X is neither of the above 
then the differential equation is 
X" + 0. 2*X' + X = 0 

The above system could be coded as 

IF(X . LT. - 0.2) THEN X' 1 + 0. 2*X' + X = 0. 2 
IF(X .GT. 0.2) THEN X M + 0.2*X' + X =-0.2 
X ,! + 0. 2*X' + X = 0. 0 
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During the execution of the solution, when the differential 
equation is referred to, the Interpreter will test the argument of each 
IF statement starting at the top. It will use the differential equation 
corresponding to the first condition it finds to be true. Note that the 
third statement is not preceded by a condition. The effect of the third 
statement is that the condition for its use is always true. Thus if the 
Interpreter does not use one of the first two differential equations it 
will use the third. 

To further clarify the selection of the right differential 
equation, consider the following: 

X" + 0. 2*X' + X = 0.0 

IF(X . LT. -0.2) THEN X" + 0. 2*X' + X =0.2 
IF(X . GT. 0.2) THEN X" + 0.2*X' + X = -0. 2 

This is the same set of statements as before except that the 
bottom statement has been moved to the top. However, the result is 
drastically different. The condition for use of the present first dif- 
ferential equation is true by default. Thus it will be used regardless 
of the value of X. Again, the point to be stressed is, only one dif- 
ferential equation can be used at a given moment; the one used will be 
the first one with a true condition. 

4. Using the Remaining Features of the Language 

Often it is convenient to be able to define and calculate some 
variables before the Differential Equation Statement. Consider the 
following coded system description: 

IF((4. 0*X + X' ) . GE. 0.0) THEN X" + 0. 02*X' + 0.3 =0.0 
IF((4. 0*X + X' ) .LT. 0.0) THEN X" + 0. 02*X' - 0.3 =0.0 

Using features of the language previously described this could be 
recoded as follows: 
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VI = + 0.3; IF((4. 0*X + X' ) . LT. 0.0) THEN VI = -VI 
X" + 0. 02 *X' + VI = 0.0 

The semicolon allows more than one statement per line. 
This recoded version has the same effect as the previous version. 
Note that there is only one differential equation. The VI used in the 
differential equation is defined in the line above in such a way that 
the effect is the same as the two differential equations of the first 
version. 

For this system an even simpler coding is possible. 

X" + 0. 02*X' + 0. 3*SIGN(4. 0*X + X') = 0.0. This example 
demonstrates the versatility of the language by showing the same 
system coded in three different ways. 
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IV. PRINCIPLES AND METHODS OF USE 



A. GENERAL COMMENTS 

The previous chapters have concentrated on various specific 
aspects of the program. In Chapter II the internal structure of the 
program was described. Chapter III described the language associated 
with use of the program. This chapter will endeavor to pull these to- 
gether with specific instructions on the use of the program. In Section 
B those areas common to both the Batch Version and the Interactive 
Graphics Version will be discussed. Sections C and D will deal with 
the peculiarities of the two versions respectively. 

B. GENERAL PRINCIPLES 

1. System Description 

One of the first tasks for the user is to encode the system 
equations in the form usable by the computer. One of the primary 
features of the program is that this is done in a very natural way. 

For simple systems this will involve writing only a single equation. 

For systems with a complicated control strategy the user may have 
to write what amounts to a small program to completely specify it. 

This general area was covered with examples in Section C of Chapter 
III. Chapter V will further amplify this with a large 
number of complete examples. 

2. Use of the Constants 

The use of the constants was described in Chapter III. How- 
ever for the completeness of this Chapter they will be mentioned 
again. Often when preparing to study a particular system it may be 
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desired to vary a parameter over several values. One way to accom- 
plish this is to vary the numerical value in the system equations. 
However this amounts to starting a whole new problem and the system 
equations must be recompiled. 

An easier way to accomplish this parameter variation is to 
use the constants. These are the constants labeled A - H. When 
writing the system equation use the constants in place of an actual 
numerical value and define them separately as desired. When it is 
desired to change the parameter, redefine the associated constants. 

Another handy use of these constants is when the same 
number is to be used several places in the system description. In 
this use a constant would be defined once and then used as often as 
desired. 



3. The Phase Plane Window 

The output of the program is a phase plane. The slopes and 
solutions are presented on this phase plane. This phase plane can be 
visualized as a window in which the user has control over the place- 
ment, size, and scales to be used in its construction. The controls 
are independent in the horizontal and vertical directions. Following 
are the control parameters and their use: 

X-SCALE The horizontal scale in units per inch. 

Y-SCALE The vertical scale in units per inch. 

X-CENTER The X coordinate of the center of the window. 

Y- CENTER The Y coordinate of the center of the window. 

X-SIZE The horizontal dimension in inches. 

Y-SIZE The vertical dimension in inches. 



40 



4. Solution Parameters 



The remaining values to be specified are the parameters of 
the solution. These include the initial time, time step, final time, 
and the initial conditions. 

In most cases, the user will let the initial time be zero and 
the final time be some positive value with an associated positive time 
step. However, this is not required. In fact the final time may be 
less than the initial time. In this later case the time step would be 
negative. 

The program can handle systems up to eighth order. Thus 
there are provisions for specifying up to eight initial conditions. 

This limitation to eighth order was arbitrarily set by the ease of 
handling the initial values on input. Eight initial values fit nicely on 
one input card. 

C. SPECIFIC DETAILS ON THE BATCH VERSION 
1 . Planning for the Batch Run 

For the batch version the user must plan the entire run 
beforehand. One may be uncertain of some parameters of the 
problem on the first run. If this is the case, the user should take 
advantage of the flexibility of this program on the first submission. 
That is, the first submission should include combinations of values 
such that when the results are returned one will know what values 
to choose for the next run. For example, take the problem of the 
choice of time step and final time. Here it would be good to make 
a guess at a suitable time step and final time. Then run an addition- 
al two solutions with these times multiplied by ten and divided by 
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ten. The added effort for these additional solutions would be trivial. 



When the results come back it will then be quite easy to decide on 
the best values. 

The same idea can apply to other areas such as scaling. The 
time required to punch the added cards for a few extra runs can save 
considerable time lost to turn-around-time if an added exploratory run 
is required. 

2. Data Card Arrangement 



Cards 1-2 


Plot title including your name (6A8) 
Format 


3 


X-SCALE(Units /inch), 

X-CENTER (Coordinate of X-CENTER) 
X-SIZE (Inches) 

Y-SCALE (As above) 

Y- CENTER (As above) 

Y-SIZE (As above) 

These six items are placed on one 
card in the above order using a 
(6F10. 0) Format. 

If (X-SIZE . GT. 9.0) then the plot is 
rotated CCW on the paper by a quarter 
turn. Due to paper size only one of 
the dimensions may exceed 9. 0 inches 
and must not be more than 15. 0 inches. 


4 


SIZE of slope markers (inches, 0.2 
works nicely), NUMBER per inch (3 
recomme nded) 

In (F10.4, 11) Format 


5, 6, 7, 8. . . 


EQUATIONS TO BE USED: Following 
the rules of Chapter III (80A1 Format). 



Then follow four cards per solution desired as 
described below: 



-First- 


VALUES of A, B, C, D, E, F, G, H (8F10. 4) 


-Second- 


Format 

Integration info (INITIAL TIME, TIME 
STEP, and FINAL TIME) (900 Solution 


- Third- 


steps Max), (3F10.4) Format. 
INITIAL CONDITIONS for X, X', X", 
. . . in (8F10. 4) Format 
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-Fourth- A "WHATS NEXT CARD" use the following 
code: 

0 - There are no requests following. 

1 - The following are the second and 

third cards with new integration 
and initial condition information 
only. I want this new solution to 
be plotted on the present graph. 

2 - The following is another four cards 

with data for another solution of the 
present system on a new graph. 

3 - The following card starts a 

completely new problem starting 
with the plot title cards. 

3. Sample Data Deck 

Following is a typical sample deck of input cards. The com- 
ments enclosed within parenthesis are not part of the cards but are 
included in this listing for clarification only. These cards accomplish 
the following. 

The system X 11 + A*X 1 + SIN(X) = 0 is compiled and two 

solutions are obtained with A = 0. 3. Then A is changed to 0. 7 and 
two more solutions are run using the same time information and 
initial conditions as with A = 0. 3. Then a completely new problem 
(a saturated servo system) is defined. 

The input cards used are: 

NELSON, H. G. Pendulum (GRAPH TITLE - 1st LINE) 



X" = 


A *X ' = SIN(X) = 0 


(GRAPH TITLE - 2nd LINE) 


2. 0 


0 . 0 


8. 0 


2. 0 0. 0 6.0 (WINDOW INFO) 


0. 2 


4 




(SLOPE MARKER SIZE AND NO. ) 


X" + A*X' + SIN(X) = 0 


(SYSTEM EQUATION) 








(END OF SYSTEM EQUATIONS) 


0. 3 






(VALUE OF A) 


0 . 0 


0. 05 


25. 0 


(TIME INFO) 


-7. 0 


3. 0 




(INITIAL CONDITIONS) 


1 






(WHAT'S NEXT SIGNAL) 


0 . 0 


0. 05 


25. 0 


(TIME INFO) 


-7. 0 


5. 0 




(INITIAL CONDITIONS) 


2 






(WHAT'S NEXT SIGNAL) 


0. 7 






(VALUE OF A) 


0 . 0 


0. 05 


25. 0 


(TIME INFO) 
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-7.0 3.0 (INITIAL CONDITIONS) 

1 (WHAT'S NEXT SIGNAL) 

0. 0 0. 05 25. 0 (TIME INFO) 

-7. 0 5. 0 (INITIAL CONDITIONS) 

3 (WHAT's NEXT SIGNAL) 

NELSON, H. G. 0902 (GRAPH TITLE - 1st LINE) 

SATURATED SERVO S YSTEM ( GRAPH TITLE - 2rd LINE) 

0.3 0.0 6.0 0.2 0.0 8.0 (WINDOW INFO) 

0. 2 4 (SLOPE MARKER SIZE AND NO. ) 

IF(X. LT. (1. 0-2. 0)*0. 2 THEN X" + 0. 2*X' + X = 0. 2 

IF(X. GT. 0. 2) THEN X" + 0. 2*X' + X = -0. 2 (SYSTEM EQUATIONS) 

X" + 0. 2*X’ + X = 0. 0 

(END OF SYSTEM EQUATIONS) 

0.0 0.05 30.0 (TIME INFO) 

0.7 (INITIAL CONDITIONS) 

(WHAT'S NEXT = BLANK = 0, THUS 
END OF DATA) 



D. SPECIFIC DETAILS ON THE INTERACTIVE GRAPHICS VERSION 

1. Planning for the Interactive Session 

Since the Interactive Version is very easy to use and, as the 
name implies, is interactive the required preparation is minimal. 

This statement assumes the first-time user has read this manual. If 
a mistake is made or a bad guess is made for a parameter value, it 
can normally be corrected before the run starts. If one notices the 
error after the run has started, there is still no problem. The pro- 
gram has a large amount of protection built into it. Hopefully, no 
mistakes or errors by the user can cause the interaction to terminate. 
Regardless of the immediate results of a mistake, control will return 
to the user for desired corrections and/or modifications. 

2. Check-off List for Loading the Program 

Following is a step by step list to be used to load and run the 
program. 

a. Turn on the XDS 9300 Computer (if not already on) using 
the instructions attached to the console. Also turn on Sense Switch 2. 
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b. Mount the binary tape of the program on one of the tape 



drives and advance it to the load point. Set the "unit select" to 0, the 
"density select" to 556 and the mode switch to automatic. The "unit 
power", "unit ready", "file protect", and "load point" lights should be 
on. 

c. Place the small card deck as specified in Appendix B in 
the card reader and press the "power on" and "start" buttons in that 
order. 

d. Ready the printer. On the XDS 9300, verify that "Sense 
Switch 2" is on. Then press the following buttons: 

HALT 

RESET 

CLEAR and CLEAR FLAGS simultaneously 

RUN 

CARDS. 

e. The card reader should begin reading cards and then the 
Tape Drive will begin operating. It will take about three minutes to 
load and link the program. 

f. Following the instructions at the AGT's, turn on the AGT 
unit to be used and load "Gated". 

g. Verify that "Gated" is ready by pressing the "Text Edit" 
button on the moveable button panel. The following should appear at 
the bottom of the screen: 

"TEXT BLOCK SELECT MODE Block 1" 

h. When the XDS console typewriter types "ENTER 
IDEV = 1 or 2" do the following: 

1. Ensure the AGT is ready 

2. Note the number of the AGT 

3. Type "IDEV = 1*", return 

or "IDEV = 2*", return as applicable. 



45 



i. This completes the preparation phase. Although this list 
appears long or complicated it is quite easy once it has been done a 
few times. Quite often the user will find both the XDS and the AGT 
already turned on and ready to go. It is strongly recommended that 
the initial use of this program be made during working hours when 
advice and/or help is available. 

3. Interacting with the Program 

The interaction starts with a question to the user. "Do you 
desire to start a completely new problem?, answer 1 YES 1 or , NO ! 
Answer "YES" or "NO" on the keyboard and then hit the return key. 

If you answer "NO", a problem already in the computer will be dis- 
played. This is useful for becoming familiar with the program. It 
is suggested for the first time user that the answer should be "NO" 
to that question. Then the new user should try various manipulations 
on the given problem until a working knowledge of the program is 
obtained. 

After each run the control returns again to this same question. 
After the first pass the answer to the question gives the user the 
option of making modifications to the present system already on the 
screen or starting a complete new problem. If there are only a few 
changes to be made it is best to selectively edit the present problem. 

If it is desired to start a new problem or make many changes, it will 
probably be faster to use the "INPUT" mode. When using the "INPUT" 
mode, the system prompts the user for the desired information. 

The easiest way to learn the system is to use it. Thus the 
following section is a step-by-step procedure for using the program 
on several specific problems. It is suggested that the user closely 
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follow this guide using the given problems. Once the user is familiar 
with the procedures for using the program and the graphics terminal 
then he can proceed on his own. 

4. Suggested First Session Example 

a. Load the program as described in Section 2. The question 
n DO YOU DESIRE TO START A COMPLETELY NEW PROBLEM? , 
ANSWER YES OR NO’ 1 will appear on the screen. 

b. Type "NO”, return. At this time the internally stored 
program will appear on the screen. In addition will be the message 
"MAKE DESIRED CHANGES AND CHOOSE NEW INITIAL CONDITIONS 
(TEXT OR LIGHT PEN) TO INITIATE NEW RUN". The signal to the 
program that the user is ready to proceed is to edit the first line of 
the initial conditions or to select initial conditions graphically with 
light pen. 



c. If satisfied with the initial conditions already on the 
screen, make the following null edit. Press the following buttons in 
succession: 

TEXT EDIT 

NEXT BLOCK - Check to see that the second line characters are 

slanted. If not use the NEXT BLOCK and/or 
PREVIOUS BLOCK until line two is slanted. The 
slanting indicates that line two is ready to be 
edited. 

GO EDIT - Corrections can now be made as desired to line 

two. For this time a null correction will be 
made, i. e. , no correction. 

TERMINATE EDIT 

The program now has all the information it needs to make a run and 
the above edit of line two has indicated that the program is ready to 
be run. When the run has terminated the question of step a. will be 
repeated. A complete loop in the program has now been completed, 

d. Again answer by typing n NO u , return. 
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e. Before this next run the goal will be to change the initial 
conditions to -3. 8 and 3. 0 respectively. To do this "select" and "go 
edit" on line two as in step c above. This time a correction will be 
made. While holding the typewriter key marked "CTRL", space over 
with the "F" (forward) key until the underscore is immediately to the 
right of the first initial condition to be corrected. Now release the 
"CTRL" key. Press the "RUB OUT" key the desired number of times 
to eliminate the portion to be changed. Then type in the new first 
initial condition. Be sure to replace as many characters as were 
removed. This keeps the whole line in the proper alignment. Next, 
space over and correct the second initial condition using the same keys 
as used to correct the first initial condition. When you are satisfied 
with the line press the terminate edit button or the "RETURN" key on 
the typewriter. Either one will terminate the edit and initiate the next 
run. The end of this run will bring you back to the question of step a. 

f. Again answer by typing "NO", return. 

g. The goal of this run is to select a set of initial conditions 
using the light pen. Press in sequence the "graphics edit", "go edit", 
and "cursor on" buttons. Insure that the light pen is turned on by 
checking for the pilot light on the box at the end of the light pen cable. 
Take the light pen, point it at the cursor, press the little white button 
on the light pen and lead the cursor around the screen. When you are 
familiar with positioning the cursor, select with the cursor a desired 
initial condition, then press the "move" and then "draw" buttons. 

With the light pen move the cursor to another initial condition on the 
screen. There should be a dot remaining at the location of the first 
initial condition pair selected. Pressing the "move", "draw" buttons 
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will select the present cursor position as a second initial condition 
pair. This can be repeated up to fifteen times if desired. When you 
are through selecting initial conditions, terminate the edit. This will, 
as when selecting the initial conditions numerically, initiate the run. 

A solution will now be obtained for each of the selected initial con- 
ditions and then control will return again to the question of step a. 

h. Again answer "NO”. 

i. This time change A to 3. 0 . To do this press "text edit" 
and then press next block until the characters in the line containing 
the A tilt. You are now ready to edit this line using the method of 
step e above. 

j. When through editing A initiate the next run using the 
technique of steps c, e, or g above. 

k. Once more answer "NO". 

l. This time it is desired to edit the system equation to 
obtain X n + X*X* + X = 0. Press the "text edit" and then the "next 
block" button about eight times until the equation characters slant 
indicating they have been selected. Make the modification as before. 
You do not need to replace as many characters as you removed when 
editing the equations. When ready initiate the next run using one of 
the previously discussed methods. 

m. When the solution has been completed and the question 
of step a returns answer "YES". In response to each question make 
responses as suggested below. After each response hit the return 
key. 
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SYSTEM EQUATIONS: 



VI = 4. 0*X + X' 

IF( V 1 . GE. 0. 0) THEN X" + 0. 02*X' + 0. 3 = 0. 0 
X" + 0. 02 *X' - 0. 3 = 0. 0 
blank line 



X-SCALE 
X- CENTER 
X-SIZE 
Y -SCALE 
Y-CENTER 
Y-SIZE 
A - H 

INITIAL TIME 
TIME STEP 
FINAL TIME 
FIRST I. C. 
remainder of I. C. 



0 . 2 
0 . 0 
6 . 0 
0 . 2 
0 . 0 
8 . 0 

hit return, will default to zero 
0 . 0 
0. 04 
10 . 0 
0 . 6 

s, hit return. 



n. Now the program will give the question "MAKE DESIRED 
CHANGES (TEXT OR LIGHT PEN) TO INITIATE NEW RUN". At 
this point you can proceed as in steps c. through 1. above as desired. 
Thus if any corrections need to be made as a result of incorrect 
responses to the above questions they can be made at this time. When 
ready to initiate the run use the method of step c. It is best to select 
the initial conditions numerically for the first run after changing the 
window parameters. Initial conditions selected with the light pen take 
on the corresponding values of that point relative to the current 
picture on the screen. 

o. You have at this time exercised all the control functions 
needed for full interactive control of the program. 
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V. APPLICATIONS TO VARIOUS PROBLEMS 



A. GENERAL COMMENTS 

The purpose of this chapter is to demonstrate the use of the 
program by using it to study several classes of control system prob- 
lems. It is not the intention of this chapter to derive the applicable 
equations, but rather to show the use of the program in their solution. 
If the purpose of the system is to reduce the error to zero or a 
minimum, then it is convenient to use error as the dependent variable 
instead of the regular output variable. 

In the following examples the problem will be defined first then 
the data cards will be shown, followed by the resultant phase plane. 
Although these examples will have been produced using the hard plots 
of the Batch Version, the input data is identical with that of the 
Interactive Graphics Version. 
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B. VARIOUS SINGLE EQUATION EXAMPLES 
Example 1 

Van der Pol's Equation 

X" - A*( 1 . 0 - X**2)*X' + B*X = 0. 0 

WHERE: 

A = 1. 0 
B = 1. 0 



PLOT PARAMETERS: 



XSCALE 


1 . 0 


YSCALE 


1 . 0 


X CENTER 


0 . 0 


Y CENTER 


0 . 0 


XSIZE 


6. 0 


YSIZE 


8. 0 



Size of slopes = 0.2 
Number of slopes per inch = 4 



THE INPUT CARDS USED: 

NELSON, H. G. 

X" - A*(l. 0 - X**2)*X' + B*X = 0. 0 

1 . 0 0.0 6.0 1.0 0.0 8.0 

0. 2 4 

X" - A*(l. 0 - X**2)*X' + B*X = 0. 0 

1.0 1.0 

0.0 0.1 20.0 

0 . 01 0 . 0 

1 

0 . 0 0 . 1 20 . 0 

1.0 0.0 

1 

0 . 0 0 . 1 20 . 0 

2.5 3.0 

1 

0 . 0 0 . 1 20 . 0 

2. 0 -3. 3 

1 

0 . 0 0 . 1 20 . 0 

-2. 5 3. 0 

1 

0. 0 0. 1 20. 0 

-2. 5 -2. 0 
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Example 1: X" - A*(l. 0 - X**2)*X' + B*X = 0. 0 




Fig. 5-1 X-SCALE = 1. 0 units /inch 
Y-SCALE = 1.0 units /inch 
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Example 2: 



X" + A*X*X' + B*X = 0. 0 

WHERE: 

A = 1. 0 
B = 1. 0 



PLOT PARAMETERS: 



XSCALE 


1. 0 


YSCALE 


1. 0 


XCENTER 


0. 0 


Y CENTER 


0. 0 


XSIZE 


6. 0 


YSIZE 


8. 0 



Size of slopes = 0.2 

Number of slopes per inch = 4 



THE INPUT CARDS USED: 

NELSON, H. G. 

X" + A*X*X' + B*X = 0. 0 
1.0 0.0 6.0 1.0 0.0 8.0 
0. 2 4 

X" + A*X*X' + B*X = 0. 0 

1.0 1.0 

0.0 0.02 15.0 

1 . 0 

1 

0.0 0.02 15.0 

2 . 0 

1 

0.0 0.02 15.0 

3. 0 

1 

0.0 0.02 15.0 

2. 5 -2. 0 

1 

0.0 0.02 15.0 

2.5 -3.5 
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Example 2: X" + A*X*X' + B*X 
A = 1. 0, B = 1. 0 



= 0 . 0 
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Example 3: 

X" +- (1. 0 - X**2)*X' + X = 0. 0 

•O 

PLOT PARAMETERS: 



XSCALE 


1. 0 


YSCALE 


1. 0 


XCENTER 


0. 0 


Y CENTER 


0. 0 


XSIZE 


6. 0 


YSIZE 


8. 0 



Size of slopes = 0.2 

Number of slopes per inch = 4 
The input cards used: 
NELSON, H. G. 



X" 


+ (1. 0 - X**2)*X' 


+ X = 0. 0 






1. 0 


0. 0 6. 0 


1. 0 


0. 0 


8. 0 


0. 2 


4 








X" 


+ (1.0- X**2)*X' 


+ X = 0. 0 







0 . 0 

0.0 0.01 8.0 

1. 0 
1 

0.0 0.01 10.0 

1. 5 
1 

0.0 0.01 8.0 

2 . 0 

1 

0 . 0 0 . 01 8 . 0 

2. 25 
1 

0.0 0.01 8.0 

-2. 05 
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Example 3: X" + (1. 0 - X**2)*X' + X = 0. 0 




Fig. 5-3 X-SCALE =1.0 units/inch 
Y-SCALE = 1. 0 units /inch 
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Example 4: 
x n + x'**2 + X = 0. 0 
PLOT PARAMETERS: 



XSCALE = 2. 0 

XCENTER = 0. 0 

XSIZE = 8. 0 

Size of slopes -= 0. 2 



Number of slopes per inch = 4 

THE INPUT CARDS USED: 
NELSON, H. G. 



X" + X 


' **2 + X = 


0 . 0 


2. 0 


0 . 0 


8. 0 


0.2 


4 




X" + X 


■ **2 + X = 


0 . 0 


0 . 0 
0 . 0 


0. 05 


40. 0 


-4. 0 


2. 5 




1 

0 . 0 


0. 05 


40. 0 


-7. 0 

1 

0 . 0 


0. 05 


40. 0 


1. 0 


3. 0 




1 

0. 0 


0. 05 


40. 0 


-5. 0 
1 

0 . 0 


0. 05 


40. 0 


-3. 0 
1 

0 . 0 


0. 05 


40. 0 


-1. o 
1 

0 . 0 


0. 05 


40. 0 


0.4 


2. 5 





YSCALE 


1. 0 


Y CENTER 


0. 0 


YSIZE 


6. 0 


0. 0 6. 0 


* 
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Example 4: X" + X'**2 + X 



o 

o 

II 




59 



Fig. 5-4 X-SCALE = 2.0 units/inch 



Example 5: 

X" + A*X + B*X**3 = 0. 0 

WHERE: 

A = -1. 0 
B = 0. 25 

PLOT PARAMETERS: 



XSCALE = 1.0 

XCENTER = 0. 0 

XSIZE = 8. 0 

Size of slopes = 0.2 



Number of slopes per inch = 4 
THE INPUT CARDS USED: 
NELSON, H. G. 



X" 


+ A*X + BrX 


— 1 


1. o 


0. 0 


8. 0 


0. 2 


4 




X" 


+ A*X + B-'-X 


*❖3 = i 


-1. o 


0. 25 




0. 0 


0. 02 


18. 0 


1. 0 
1 


1.2 




0. 0 


0. 02 


18. 0 


1. o 
1 


2. 0 




0. 0 


0. 04 


30. 0 


0. 0 

1 


0. 03 




0. 0 
0. 3 

1 


0. 02 


18. 0 


0. 0 
-0. 3 
1 


0. 02 


18. 0 


0. 0 
1. o 
1 


0. 02 


18. 0 


0. 0 
-1.0 
1 


0. 02 


18. 0 


0. 0 
-1. 5 


0. 02 


18. 0 


i. o 

1. 5 


0. 02 


18. 0 



YSCALE = 


1 . 0 


Y CENTER 


0 . 0 


YSIZE 


6. 0 



0 . 0 6 . 0 
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Example 5: X M + A*X + B*X**3 
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Fig. 5-5 X-SCALE =1.0 units/inch 
Y-SCALE = 1. 0 units/inch 



Example 6: 

X" + (1. 0 - ABS(X))*X* + X = 0. 0 
PLOT PARAMETERS: 



XSCALE 


1 . 0 


YSCALE 


— 


1 . 0 


XCENTER 


0 . 0 


Y CENTER 




0 . 0 


XSIZE 


6. 0 


YSIZE 


— 


8. 0 



Size of slopes = 0. 2 
Number of slopes per inch = 4 



THE INPUT CARDS USED: 

NELSON, H. G. 

X" + (1. 0 - ABS(X))*X' + X = 0. 0 

1.0 0.0 6.0 1.0 0.0 8.0 

0. 2 4 

X" + (1.0- ABS(X))*X' + X = 0. 0 



0 . 0 

0 . 0 


0. 02 


15. 0 


1 . 0 
1 

0 . 0 


0. 02 


15. 0 


1. 5 
1 

0 . 0 


0. 02 


15. 0 


2. 0 

1 

0 . 0 


0. 02 


15. 0 


2. 8 
1 

0 . 0 


0. 02 


15. 0 


-2. 5 
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Example 6: X" + (1. 0 - ABS(X))*X' + X = 0. 0 




Fig. 5-6 X-SCALE = 1. 0 units/inch 
Y-SCALE =1.0 units /inch 
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Example 7: 



X" + X' + X*ABS(X) = 0. 0 
PLOT PARAMETERS: 



XSCALE 


2, 0 


YSCALE 


= 2. 0 


XCENTER 


0 . 0 


Y CENTER 


= 0 . 0 


XSIZE 


6. 0 


YSIZE 


= 8. 0 



Size of slopes = 0. 2 
Number of slopes per inch = 4 



THE INPUT CARDS USED: 

NELSON, H. G. 

X" + X' + X*ABS(X) = 0. 0 

2 . 0 0 . 0 6 . 0 2 . 0 0 . 0 8 . 0 

0. 2 4 

X" + X' + X*ABS(X) = 0. 0 



0. 0 






0. 0 
3. 8 
1 


0. 01 


8. 0 


0. 0 
-3. 0 
1 


0. 01 


8. 0 


0. 0 
5. 0 
1 


0. 01 


8. 0 


0. 0 
-5. 5 


0. 01 


8. 0 
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Example 7: X" + X' + X*ABS(X) = 0. 0 




Fig. 5-7 X-SCALE = 2. 0 units/inch 
Y-SCALE = 2. 0 units /inch 
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C. VARIOUS TYPES OF PENDULUM EXAMPLES 
Example 8: 

Linear approximation to the pendulum, g/1 = 1 
No Damping 

X" + X = 0 



PLOT PARAMETERS: 



XSCALE 


2. 0 


YSCALE 


2. 0 


XCENTER 


0 . 0 


Y CENT ER 


0 . 0 


XSIZE 


6. 0 


YSIZE 


8. 0 



Size of slopes = 0. 2 
Number of slopes per inch = 4 



THE INPUT CARDS USED: 
NELSON, H. G. PENDULUM 



X" 


+ X = 0 




2. 0 
0. 2 
X" 


0. 0 
4 

+ X = 0 


6. 0 


0. 0 

2. 0 

1 


0. 02 


10. 0 


0. 0 

1. 0 

1 


0. 02 


10. 0 


0. 0 
4. 0 

1 


0. 02 


10. 0 


0. 0 
5. 5 


0. 02 


10. 0 
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Example 8: X" + X = 0 




Y-SCALE = 2. 0 units /inch 
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Example 9: 



Pendulum without damping, g/1 = 1 
X" + SIN(X) = 0 
PLOT PARAMETERS: 



XSCALE 


2. 0 


YSCALE 


2. 0 


XCENTER 


0 . 0 


Y CENTER 


0 . 0 


XSIZE 


8. 0 


YSIZE 


6. 0 



Size of slopes = 0. 2 
Number of slopes per inch = 4 
THE INPUT CARDS USED: 



NELSON, H. G. 


PENDULUM 


X" 


+ SIN(X) = 0 






2. 0 


0. 0 


8. 0 


2. 0 


0. 2 


4 






X" 


+ SIN(X) = 0 






0. 0 


0. 02 


10. 0 




-6. 0 


1. 0 






1 








0. 0 


0. 02 


15. 0 




-6. 0 


2. 0 






1 








0. 0 


0. 02 


10. 0 




1 

G N 

O 


3. 0 






l 








0. 0 


0. 02 


10. 0 




1 

O 


4. 0 






l 








0. 0 


0. 02 


10. 0 




0. 0 


1. 0 






1 








0. 0 


0. 02 


' 10. 0 




6. 0 


1. 0 






1 








0. 0 


0. 02 


15. 0 




6. 0 


-2. 0 






1 








0. 0 


0. 02 


10. 0 




6. 0 


-3. 0 






1 




. 




0. 0 


0. 02 


10. 0 




6. 0 


-4. 0 
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Example 9: X" + SIN(X) 
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Fig. 5-9 X-SCALE = 2. 0 units/inch 



Example 10: 



Pendulum with velocity damping. 
X" + A*X'+ SIN(X) = 0 
WHERE: 

A = 0. 3 



PLOT PARAMETERS: 



XSCALE 


2. 0 


YSCALE 


2. 0 


XCENTER 


0 . 0 


Y CENTER 


0 . 0 


XSIZE 


8. 0 


YSIZE 


6. 0 



Size of slopes = 0. 1 
Number of slopes per inch = 4 



THE INPUT CARDS USED: 

NELSON, H. G. PENDULUM 
X" + A*X' + SIN(X) = 0 

2 . 0 0 . 0 8 . 0 2 . 0 0 . 0 6 . 0 

0. 1 4 

X" + A*X' + SIN(X) = 0 



0. 3 
0. 0 


0. 05 


25. 0 


-7. 0 


3. 0 




1 

0. 0 


0. 05 


25. 0 


-7. 0 


5. 0 




1 

0. 0 


0. 05 


25. 0 


-7. 0 


2. 0 




1 

0. 0 


0. 05 


25. 0 


6. 0 


-5. 0 
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Example 10: X" +A*X' + SIN(X) 
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Fig. 5-10 X-SCALE = 2. 0 units/inch 



Example 11: 

Pendulum with velocity squared damping. 
X" + A*X'*ABS(X') + SIN(X) = 0 
WHERE: 

A = 1. 0 

PLOT PARAMETERS: 



XSCALE 


2. 0 


YSCALE 


2. 0 


X CENTER 


0 . 0 


Y CENTER 


0 . 0 


XSIZE 


8. 0 


YSIZE 


6. 0 



Size of slopes = 0. 2 
Number of slopes per inch = 4 



THE INPUT CARDS USED: 

NELSON, H. G. PENDULUM 
X" + A*X'*ABS(X') + SIN(X) = 0 
2 . 0 0 . 0 8 . 0 2 . 0 0 . 0 6.0 
0. 2 4 

X" + A*X'*ABS(X') + SIN(X) = 0 



1. 0 
0. 0 


0. 03 


20. 0 


1. 0 


4. 0 




1 

0. 0 


0. 03 


20. 0 


-4. 0 


4. 0 




1 

0. 0 


0. 03 


2 0. 0 


-7. 0 


4. 0 




1 

0. 0 


0. 03 


20. 0 


5. 0 


-4. 0 




1 

0. 0 


0. 03 


20. 0 


2. 0 


4. 0 




1 

0. 0 


0. 03 


2 0. 0 


-2. 5 


-4. 0 
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Example 11: X" + A*X'*ABS(X') + SIN(X) 
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Fig. 5 -11 X-SCAL.E = 2. 0 units/inch 



D. SATURATED SYSTEMS 



Example 12: 

A saturated servo system. Saturation occurs when absolute error 
signal exceeds 0. 2. 

IF(X. LT. -0. 2) THEN X" + 0. 2*X' + X = 0. 2 
IF(X. GT. 0. 2) THEN X" + 0. 2*X' + X = -0. 2 
X" + 0. 2*X' + X = 0. 0 



PLOT PARAMETERS: 



XSCALE 


0. 3 


yscale 


= 


0. 2 


XCENTER 


0 . 0 


Y CENTER 


= 


0 . 0 


XSIZE 


6. 0 


YSIZE 


— 


8. 0 



Size of slopes = 0. 2 
Number of slopes per inch = 4 



THE INPUT CARDS USED: 

NELSON, H. G. 0902 

SATURATED SERVO SYSTEM 
0 . 3 0 . 0 6 . 0 0 . 2 0 . 0 8 . 0 
0. 2 4 

IF(X. LT. (1. 0-2. 0)*0. 2) THEN X" + 0. 2*X' + X = 0. 2 
IF(X. GT. 0. 2) THEN X" +0. 2*X' + X = -0. 2 
X" + 0. 2*X' + X = 0. 0 



0.0 0.05 30.0 

0. 7 
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Example 12: Saturated Servo System 




Fig. 5-12 X-SCALE = 0. 3 units/inch 
Y-SCALE =0.2 units /inch 
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E. RELAY SYSTEMS 



Example 13: 

IDEAL RELAY 

X" + 0. 3*X' + 1. 5*SIGN(X) = 0. 0 
PLOT PARAMETERS: 



XSCALE 


1. 0 


YSCALE 


= 0. 6 


XCENTER 


0. 0 


YCENTER 


II 

O 

o 


XSIZE 


6. 0 


YSIZE 


= 8. 0 



Size of slopes = 0.2 
Number of slopes per inch = 4 



THE INPUT CARDS USED: 

NELSON, H. G. 

IDEAL RELAY 

1 . 0 0.0 6.0 0.6 0.0 8.0 

0.2 4 

X" + 0. 3*X' + 1. 5*SIGN(X) = 0. 0 



0. 0 0. 1 85. 0 

2 . 8 
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Example 13: X" + 0. 3*X' + 1. 5*SIGN(X) = 0. 0 




Fig. 5-13 X-SCALE - 1. 0 units /inch 
Y-SCALE = 0. 6 units /inch 
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Example 14: 



Ideal Relay with Rotated Switching Line 
X" + 0. 2*X' + 1. 0*SIGN(4. 0*X - X ! ) = 0. 0 
PLOT PARAMETERS: 



XSCALE 


2. 0 


YSCALE 


1. o 


XCENTER 


0. 0 


Y CENTER 


0. 0 


XSIZE 


6. 0 


YSIZE 


8. 0 



Size of slopes = 0. 2 
Number of slopes per inch = 4 



THE INPUT CARDS USED: 

NELSON, H. G. 

IDEAL RELAY WITH ROTATED SWITCHING LINE 



2. 0 


0. 0 6. 0 1.0 


O 

o 


0. 2 


4 




X" 


1 + 0. 2*X' + 1. 0*SIGN(4. 0*X 


- X') = 0 


0. 0 


0. 0 
5. 5 


0. 1 80. 0 
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Example 14: X" + 0. 2*X' + 1. 0*SIGN(4. 0*X - X') = 0. 0 




Y-SCALE = 1. 0 units /inch 
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Example 15: 



Ideal Relay with Rotated Switching Line 
X" +0. 2*X' + 1. 0*SIGN(4. 0*X + X') = 0. 0 
PLOT PARAMETERS: 



XSCALE 


2. 0 


YSCALE 


00 

o 


XCENTER = 


o 

» 

o 


Y CENTER 


o 

o 


XSIZE 


6. o 


YSIZE 


o 

00 



Size of slopes = 0. 2 
Number of slopes per inch = 4 



THE INPUT CARDS USED: 

NELSON, H. G. 

IDEAL RELAY WITH ROTATED SWITCHING LINE 

2.0 0.0 6.0 0.8 0.0 8.0 

0. 2 4 

X" +0. 2*X' + 1. 0*SIGN(4. 0*X + X' ) = 0. 0 

o 

0 . 0 

0 . 0 0 . 1 20 . 0 

5. 5 
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Example 15: X" + 0. 2*X' + 1. 0*SIGN(4. 0*X + X') = 0. 0 




Fig. 5-15 X-SCALE = 2. 0 units/inch 
Y-SCALE = 0. 8 units /inch 
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Example 16: 



Ideal Relay with Rotated Switching Line 
X" + 0. 2*X' + 1. 0*SIGN(2. 0 *X + X') = 0. 0 
PLOT PARAMETERS: 



XSCALE 


2. 0 


YSCALE 


1 . 0 


X CENTER 


0 . 0 


YCENTER = 


0 . 0 


XSIZE 


6. 0 


YSIZE 


8. 0 



Size of slopes = 0. 2 
Number of slopes per inch = 4 



THE INPUT CARDS USED: 

NELSON, H. G. 

IDEAL RELAY WITH ROTATED SWITCHING LINE 

2.0 0.0 6.0 1.0 0.0 8.0 

0. 2 4 

X" + 0. 2*X' + 1. 0*SIGN(2. 0*X + X') = 0. 0 

0 . 0 

0 . 0 0 . 1 20 . 0 

5. 5 
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Example 16: X" + 0. 2*X' + 1. 0*SIGN(2. 0*X + X') = 0. 0 




Fig. 5-16 X-SCALE = 2. 0 units /inch 
Y-SCALE = 1. 0 units /inch 
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Example 17: 



Ideal. Relay with Rotated Switching Line 
X" + 0.2*X' + 1. 0*SIGN(2. 0*X - X') = 0. 0 



PLOT PARAMETERS: 



XSCALE 


2. 0 


YSCALE 


= 0. 8 


XCENTER 


0 . 0 


Y CENTER 


= 0 . 0 


XSIZE 


6. 0 


YSIZE 


= 8. 0 



Size of slopes = 0. 2 
Number of slopes per inch = 4 



THE INPUT CARDS USED: 

NELSON, H. G. 

IDEAL RELAY WITH ROTATED SWITCHING LINE 

2.0 0.0 6.0 0.8 0.0 8.0 

0. 2 4 

X" + 0. 2*X' + 1. 0*SIGN(2. 0*X - X') = 0. 0 

0 . 0 

0. 0 0. 1 80. 0 

5. 5 
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Example 17: X" + 0. 2*X' + 1. 0*SIGN(2. 0*X - X') = 0.0 




Fig. 5-17 X-SCALE = 2. 0 units /inch 



Y-SCALE = 0. 8 units /inch 
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F. DEAD ZONE 



Example 18: 

IF(X. LT. -A. OR. X. GT. A) THEN X" + B*X' + C*X = C*A*SIGN(X) 
x ti + B*X' = 0. 0 

WHERE: 

A = 0. 3 
B = 0. 2 
C = 1. 0 

PLOT PARAMETERS: 



XSCALE 


0. 5 


YSCALE 


Z2 


0. 3 


XCENTER 


0. 0 


Y CENTER 


- 


0. 0 


XSIZE 


6. 0 


YSIZE 


- 


* 8. 0 



Size of slopes = 0. 2 
Number of slopes per inch = 4 



THE INPUT CARDS USED: * 

NELSON, H. G. 

DEAD ZONE 

0. 5 0. 0 6. 0 0. 3 0. 0 8. 0 

0. 2 4 

IF(X. LT. -A. OR. X. GT. A) THEN X" + B*X' + C*X = C*A*SIGN(X) 



X" 


+ B*X' 


= 0. 0 


0. 3 


0. 2 


1. 0 


0. 0 


0. 1 


50. 0 


0. 0 

1 


1. 3 




0. 0 
1. 0 


0. 1 


50. 0 
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Example 18: Dead Zone 

IF(X. LT. -A. OR. X. GT. A) THEN X" + B*X' + C*X = C*A*SIGN(X) 
X" + B-X' = 0. 0 
A = 0. 3, B = 0. 2, C = 1. 0 




Fig. 5-18 X-SCALE = 0. 5 units/inch 
Y-SCALE = 0. 3 units/inch 
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Example 19: Dead Zone 

IF(ABS(X). LT. 2. 0) THEN X" 
X" +-0. 2*X' + X = 0. 0 

PLOT PARAMETERS: 



XSCALE = 2. 0 

XCENTER = 0. 0 

XSIZE = 8. 0 

Size of slopes = 0. 2 



Number of slopes per inch = 4 



+ 0. 2*X' = 0. 0 



YSCALE 


— 


2. 0 


YCENTER 


= 


0 . 0 


YSIZE 


- 


6. 0 



THE INPUT CARDS USED: 

NELSON, H. G. 

DEAD ZONE 

2 . 0 0.0 8.0 2.0 0.0 6.0 

0. 2 4 

IF(ABS(X). LT. 2. 0) THEN X" + 0. 2*X' = 0. 0 



X" 


+ 0. 2*X' 


+ X = 0. 0 


0 . 0 


0 . 0 
6. 0 


0. 05 


35. 0 



88 



Example 19: IF(ABS(X). LT. 2. 0) THEN X" + 0. 2*X 
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G. A SECOND ORDER LINEAR SYSTEM 



Example 20: 

The following three plots demonstrating the effect of various 
damping were the result of the input cards listed below. 

X" +2. 0*A*B*X' + B**2*X = 0. 0 
Where: 

A = is the natural frequency 
B = is the damping factor 

PLOT PARAMETERS: 



XSCALE 


2. 0 


YSCALE 


= 2. 0 


XCENTER 


0 . 0 


Y CENTER 


= 0 . 0 


XSIZE 


8. 0 


YSIZE 


= 6. 0 



Size of slopes = 0.2 
Number of slopes per inch = 4 



THE INPUT CARDS USED: 

NELSON, H. G. 

X" +2. 0*A*B*X' + B**2*X = 0. 0 
2 . 0 0 . 0 8 . 0 2 . 0 0 . 0 6 . 0 
0. 2 4 

X" +2. 0*A*B*X' + 8**2*X = 0. 0 

1.0 1.0 

0.0 0.1 10.0 

7. 0 
1 

0.0 0.1 10.0 

-7. 0 
1 

0. 0 0. 1 15. 0 

3. 5 
1 

0.0 0.1 10.0 

-3. 5 
2 

0.5 1.0 

0. 0 0. 1 15. 0 

7. 0 
1 
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0 . 0 
- 7 . 0 
1 


0 . 1 


15 . 0 


0 . 0 
3 . 5 
1 


0 . 1 


15 . 0 


0 . 0 
- 3 . 5 
2 


0 . 1 


15 . 0 


0 . 1 


1 . 0 




0 . 0 
5 . 0 


0 . 1 


80 . 0 
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Example 20(a) X" + 2. 0*A*B*X' + B**2*X 
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Fig. 5-20(a) X -SCALE = 2. 0 units/inch 



Example 20(b) X" + 2. 0*A*B*X' + B**2*X 
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Fig. 5-20(b) X -SCALE = 2. 0 units/inch 



Example 20(c) X" + 2. 0*A*B*X' + B**2*X 



o 

o 

II 




94 



Fig. 5 -20(c) X-SCAL.E = 2. 0 units/inch 
Y-SCALE = 2. 0 units /inch 



H. A FIFTH ORDER SYSTEM 



Example 21: 

IF(T .GT. 0.0015) SKIP 3 
VI = 3000511. 6 + A 
V2 = 285120. 0 + A *35 1 1 . 0 
V3 = A*38500. 0 

X'"" + 122. 2 *X 1 " 1 + 50835. 8*X"' + V1*X" + V2*X' + V3*X = 0. 0 



WHERE: 

G = 30212. 0 

PLOT PARAMETERS: 



XSCALE 




2. 0 


YSCALE 


— 


20. 0 


XCENTER 


= 


0 . 0 


Y CENTER 


- 


0 . 0 


XSIZE 


- 


6. 0 


YSIZE 


= 


• 

8. 0 



THE INPUT CARDS USED: 

NELSON, H. G. 

FIFTH ORDER SYSTEM , 

2 . 0 0 . 0 6 . 0 20 . 0 0 . 0 8 . 0 

IF(T. GT. 0. 0015) SKIP 3 
VI = 3000511. 6 + A 
V2 = 285120. 0 + A *3 511. 0 
V3 = A *38500. 0 

X"'" + 122. 2*X"" + 50835. 8*X'" + V*X" + V2*X' + V3*X = 0. 0 



30212. 0 
0. 0 


0. 001 


0. 4 


-3. 0 
1 

0. 0 


0. 001 


0. 4 


2. 5 

1 

0. 0 


0. 001 


0. 4 


-5. 1 
1 

0. 0 


0. 001 


0. 4 


4. 6 
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Example 21: A FIFTH ORDER SYSTEM 

IF(T . GT. 0.0015) SKIP 3 
VI = 3000511. 6 + A 
V2 ^285120. 0 + A *3511. 0 
V3 = A *38500. 0 

X'"" + 122. 2 *X 1 1 " + 50835. 8*X ,M + V1*X" + V2*X' + V3*X = 0. 0 




Fig. 5-21 X-SCALE = 2. 0 units /inch 
Y-SCALE = 20. 0 units /inch 
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I. HYSTERESIS AND DEAD ZONE 



Example 22: 

RELAY WITH DEAD ZONE AND HYSTERESIS 

THE SYSTEM DESCRIPTION EQUATIONS AND/OR PROGRAM: 

V3 = T 

IF(ABS(X) . LE. A)V1 = 0. 0 
IF(ABS(X) . GT. B)V1 = SIGN(X) 

X" + 0. 2*X' + VI = 0. 0 
WHERE: 

A = 0. 8 
B = 1. 2 

The first statement is used to suppress the slopes. 

PLOT PARAMETERS: 



XSCALE 


= 


1 . 0 


YSCALE 


= 0. 5 


XCENTER 


= 


0 . 0 


YCENTER 


= 0 . 0 


XSIZE 


- 


6. 0 


YSIZE 


= 8. 0 



THE INPUT CARDS USED: 

NELSON, H. G. 

RELAY WITH DEAD ZONE AND H YSTERISES 
1.0 0. 0 6. 0 0. 5 0. 0 8. 0 

V3 = T 

IF(ABS(X). LE. A) VI = 0. 0 
IF(ABS(X). GT. B) VI = SIGN(X) 

X" + 0. 2*X' + VI = 0. 0 

0.8 1.2 

0. 0 0. 04 30. 0 

- 2 . 8 



97 



Example 22: RELAY WITH DEAD ZONE AND HYSTERESIS 




Fig. 5-22 X-SCALE 



1. 0 units / inch 



Y-SCALE 



0. 5 units /inch 
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J. NON-LINEAR RESTORING FORCES 
Example 23: 

SECOND ORDER SYSTEM WITH A RESTORING FORCE CORRESPONDING 
TO A NON-LINEAR SPRING 

THE SYSTEM DESCRIPTION EQUATION: 

X" + 0. 9*X' + (X + A*X**3) = 0. 0 

WHERE: 

A = 0. 2 



PLOT PARAMETERS: 



XSCALE 


2. 0 


YSCALE 


2, 0 


XCENTER 


0 . 0 


Y CENTER 


0 . 0 


XSIZE 


6. 0 


YSIZE 


8. 0 



Size of slopes = 0. 2 
Number of slopes per inch = 4 
THE INPUT CARDS USED: 

NELSON, H. G. 

SECOND ORDER SYSTEM WITH NON-LINEAR SPRING 



2. 0 


0. 0 


6. 0 


0. 2 


4 




X" + 


0, 9*X' 


+ (X + A*X 


0. 2 
0. 0 


0. 02 


15. 0 


-4. 0 
2 

-0. 1 
0. 0 


0. 02 


15. 0 


-3. 0 
1 

0. 0 


0. 02 


15. 0 


1. o 


-7. 0 




1 

0. 0 


0. 02 


15. 0 


-5. 5 


5. 0 




1 

0. 0 


0, 02 


15. 0 


-5. 2 


5. 0 




1 

0. 0 


0. 02 


15. 0 


5. 0 


-5. 0 
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Example 23(a): X" + 0. 9*X' + (X + A*X**3) = 0. 0 



A = -0. 1 (gives a "soft" spring) 




Y-SCALE =2.0 units /inch 
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Example 23(b): X" + 0. 9*X' + (X + A*X**3) = 0. 0 




Y-SCALE = 2. 0 units /inch 
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Example 24: 

A FORCED SYSTEM WITH A NON-LINEAR RESTORING FORCE 

THE SYSTEM DESCRIPTION: 

X" + 0. 9*X' + (X + 0. 2*X**3) = A*COS(B*T) 

WHERE: 

A = 3. 0 
B = 0. 7 

PLOT PARAMETERS: 



XSCALE 


1. 0 


YSCALE 


II 

►— * 
o 


XCENTER = 


0. 0 


YCENTER 


= 0. 0 


XSIZE 


6. 0 


YSIZE 


= 8. 0 



THE INPUT CARDS USED: 

NELSON, H. G. 

FORCED NON-LINEAR SECOND ORDER SYSTEM 

1.0 0,0 6.0 1.0 0.0 8.0 

X" + 0. 9*X' + (X + 0. 2*X**3) = A*COS(B*T) 

3.0 0. 7 

0.0 0.04 30.0 

0.0 0.0 

3 

NELSON, H. G. 

FORCED NON-LINEAR SECOND ORDER SYSTEM 



2. 0 


0 . 0 


6. 0 2.0 0 


X" + 


(X + 0. 


2*X**3) = A*COS(B*T) 


5. 0 


0. 8 




0 . 0 


0. 04 


30. 0 


0 . 0 


0 . 0 




2 






5. 0 


1 . 0 




0 . 0 


0. 04 


30. 0 


0 . 0 


0 . 0 
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Example 24(a): X" + 0. 9*X' + (X + 0. 2*X**3) = A*COS(B*T) 

A = 3. 0, B = 0. 7 




Fig. 5-24(a) X-SCALE = 1 . 0 units /inch 
Y-SCALE = 1.0 units/inch 
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Example 24(b): X" + 0. 0*X' + (X + 0. 2*X**3) = A*COS(B*T) 




Fig. 5-24(b) X-SCALE = 1.0 units/inch 
Y-SCALE = 1. 0 units /inch 
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Example 24(c): X" + 0. 0*X' + (X + 0. 2*X**3) = A*COS(B*T) 




Fig. 5-24(c) X-SCALE = 1. 0 units /inch 
Y-SCALE = 1. 0 units /inch 
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K. BANG-BANG SYSTEMS 



Example 25. 

A SECOND ORDER BANG-BANG SYSTEM 
SYSTEM EQUATION: 

X" + SIGN(X + 0. 5*X'*ABS(X')) = 0. 0 
PLOT PARAMETERS: 



XSCALE 


2. 0 


YSCALE 


1. 0 


XCENTER 


0. 0 


Y CENTER = 


0. 0 


XSIZE 


O 

• 

00 


YSIZE 


6. 0 



Size of slopes = 0. 2 
Number of slopes per inch = 4 



THE INPUT CARDS USED: 

NELSON, H. G. 

SECOND ORDER BANG-BANG 

2.0 0.0 8.0 1.0 0.0 6.0 

0. 2 4 

X" + SIGN(X + 0. 5*X'*ABS(X')) = 0. 0 



0.0 0.02 15.0 

-7. 0 
1 

0.0 0.02 15.0 

7. 0 
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Example 25: A SECOND ORDER BANG-BANG SYSTEM 




n 

W 

< 

u 

co 

i 



107 



, 0 units / inch 



Example 26. 



A THIRD ORDER BANG- BANG SYSTEM 
SYSTEM EQUATIONS: 

VI = SIGN(X* + 0. 5*X"*ABS(X")) 

V2 = X + (1. /3. )*X"**3 + Vl*X"*X' 

V3 = VI *(0. 5*X"**2 + V1*X')**1. 5 
X'" + SIGN(V2 + V3) = 0. 0 
PLOT PARAMETERS: 



XSCALE 


2. 0 


YSCALE 


1. o 


XCENTER 


0. 0 


Y CENTER 


0. 0 


XSIZE 


6. 0 


YSIZE 


8. 0 



The input cards used: 

NELSON, H. G. 

THIRD ORDER BANG-BANG SYSTEM 

2 . 0 0 . 0 6 . 0 1 . 0 0 . 0 8 . 0 

VI = SIGN(X' + 0. 5*X"*ABS(X")) 

V2 = X + (1. /3. )*X"**3 + V1*X"*X 
V3 = V1*(0. 5*X"**2 + V1*X')**1. 5 
X"' + SIGN(V2 + V3) = 0. 0 
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Example 26: A THIRD ORDER BANG-BANG SYSTEM 




Fig. 5-26 X-SCALE = 2. 0 units /inch 
Y-SCALE = 1. 0 units /inch 
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VI. CONCLUSIONS AND SUGGESTIONS FOR FURTHER RESEARCH 



This thesis has introduced an easy to use, flexible, and powerful 
analysis tool for the study of non-linear ordinary differential equations. 
The primary application emphasis has been in relation to control 
system problems for which many examples have been presented. It 
has been shown that the use of an automatic built-in compiler coupled 
with interactive graphics result in a very useful analysis program. 

For a large number of quite diverse examples it has been demonstrated 
in this thesis that one can very adequately specify a problem using 
data cards or by typing on an interactive graphics screen. 

This study has suggested several basic ideas for further research. 

1. An extension of the program following the development presented 
in Chapter II would give the program the capability of solving systems 
described by several simultaneous non-linear differential equations. 

2. If one were concerned with execution efficiency, it would be 
very worthwhile to rewrite the interpreter in assembly language. The 
interpreter is called four times for each solution point and once for 
each slope line. Thus the slope generation and solution generation 
times could be speeded up considerably with a faster interpreter* The 
other parts of the program are used relatively infrequently and there- 
fore efficiency is not an important factor. 

3. It would be worthwhile to provide a choice of integration 
schemes to the user in addition to the present fourth order Runge- 
Kutta. 



no 



4. A built-in compiler could be used in several other applications. 
One of the very obvious areas is that of optimization. For example, 

a very useful directed search algorithm called "DIRECT" is available 
on the XDS 9300. The function to be minimized is incorporated into a 
function subroutine which must be supplied by the user. A modification 
of the compiler and interpreter presented in this thesis could be used 
as a very powerful interface for describing the function and its 
boundaries. 

5. It would be very worthwhile developing methods to assist in 
the analysis of systems higher than second order. In general, the 
presentation of the grid of slopes is not worthwhile for systems of 
higher than second order since slope lines are in a space which has a 
dimension equal to the order of the differential equation. 
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APPENDIX A 



REQUIRED JOB CONTROL LANGUAGE FOR BATCH VERSION 



1. When using the compiled binary object deck use the following 
procedure: 



// JOB CARD -Green- 
// EXEC FORTLG, REGION. GO=160K 
//LINK. SYSIN DD * 



BINARY OBJECT DECK 

//GO. FT06F001 DD S YSOUT=A, SPACE=( CYL, (6)) 

/ / GO. S YS PLOT DD SYSOUT = C, SPACE=(TRK, (1, 10)) 
//GO. SYSPLOTS DD UNIT=S YSDA, SPACE=(TRK, (2, 6)) 
//GO. SYSIN DD * 



DATA 



/# -Orange- 



2. If using the source deck use the following procedure: 

/ / JOB CARD -Green- 

// EXEC FORTCLGP, REGION. FORT = 150K, REGION. GO=160K 
//FORT. SYS PRINT DD SYSOUT=A, SPACE=(CYL, (1, 1)) 

//FORT. SYSIN DD * 

FORTRAN SOURCE PROGRAM 



/* 

//GO. FT06F001 DD SYSOUT=A, SPACE(CYL, (6)) 
//GO.SYSPLOT DD SYSOUT = C, SPACE=(TRK, (1, 10)) 
//GO. SYSPLOTS DD UNIT=S YSDA, SPACE=(TRK, (2, 6)) 
//GO. SYSIN DD * 

DATA 



/* - Orange - 
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APPENDIX B 



1. Required control cards for Interactive Graphics Version. 

BOOT CARD 

PATCH DECK 

6 AGT 
6 JOB 

A ASSIGN BI =MTOA 
A ASSIGN GO =MTOA 
A LOAD XM, MAP 

A SEG MAIN-(SINIT, SSETUP2, SCOMP, (SGRID-SINTER ), (SSOLVE- 
ADATA SINTER)) 
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THE DATA CARDS SHOULD BE ARRANGEO AS FOLLOWS: 



PROGRAM A 



LISTING OF BATCH VERSION 
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24 CALL SOLVE(TI,DEL»TF,XI»YI,ITITLE) 
IP(LSTAT.bQo3 ) GO TO 31 
RcAO ( !> » 10 ) NEXT 
10 FORMAT (II) 

CALL SETIMF 
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SUBROUTINE SNORM ( SCALE ) 
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V ( J ) = - Y ( J ) 

CALL DRAWt 1C* Y, X, 1 , n * LABEL , I TITLE, YSCALE , XSCALE, IYRIGH* I XUP,2*2 



❖ NXSIZE f NYS IZE , 1 , LAST ) 
RETURN 
END 
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92 Y ( J ) = -Y(J) 

CALL DR AW (2 ,Y, X»2,G, LAB EL ( L) , ITITLEtC. ,C. f 0. ,M t 0,0, 0,0,0, LAST) 

1 CONTINUE 
99 RETURN 

END 
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I F ( ?oC- I >15,15,16 
WRITE (6, 10 )T 

FORMAT ( • SOLUTION TERMINATED AT TIME =',FK.5 
' WITH 900 SOLUTION POINTS') 
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FORTRAN 4 VERSION OF RUNGF-KUTT A- G I LL ROUTINE 
X t Y t F i H ARE DOUBLE-PRECISION. MAX N = 25 
FUNCTION RKLDEQ ( N » Y , F , X , H, NT ) 

DOUBLE PRECISION Y , F , X , H , Q , HI , H2 , H3 , H6 
DIMENSION Y ( 2 ) » F< ?) , Q( 25) 
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DOUBLE PRECISION 7 T I K ) , ZDOT(K), TT , ZL ARGE 
COP NON ZT, /DOT, TT, L, N, LFUNCT 
CONMON/COL F/A ,B,C,D,E,F,G,H 
COMMON / S T ATUS/I. STAT 
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170 WRITE ( 6,2220) 

2220 FORMAT ( • UNBALANCED PARENS IN LOGICAL EXPRESSION, PLEASE CORRECT' 
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CHECK FOR EQUAL SIGN AND TAKE APRCPRIATE ACTION TF FOUND 
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IF( INPUT( P 11 + 1 ) „iME.LO) GO TO 162 
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RECOGNIZING NUMBERS 
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WRITE ( 6 » I 3' ' 1 ) 

1801 FORMAT ( * ^FOL LOWING IS THE CODED FQUIVALENT OF THE INPUT STRING: 
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GO TO 191 

187 IF( INPUT (Pll) • NE • N7 ) GO T9 188 
INTEG = 7 

G9 TO 191 

188 IF( INPUT(PII) • NE • N8) GO TO 189 
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195 IF(STAC’<A(PSAT-1)*NE.1E5) 30 T0 194 
IF (STACKA( PSAT ).NE«122) G0 T0 196 
STACKA(PSAT-1 ) = 127 
PSAT = PSAT - 1 
GO T 0 194 
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I F ( P I 1 »LT • 96 ) GO TO 761 

WRITE ( 3 , 1 141 ) 

1141 FORMAT (100X, » FIRST PASS 0<’) 

GO TO SOI 

761 IF( IRPUT (PI 1 ) .EQ. BLANK) GO TO 157 
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WRITE (3/ 1800) ( STACKA ( I ) t I » PSAT ) 
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ENCODE ( 96 / 11 # 1 TF XT ( 1 # 30 ) ) 

CALL TEXTS ( I DEV# I TEXT (1 / 30 ) / ?4, 30# 1 # 1 , I ER ) 

RETURN 

END 
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THIS SUBROUTINE SOLVES THE ACTUAL SYSTEM AND SENDS THE S8LUT I 0N 
GUT T6 THE ADAGE AND ALS8 GIVES HARD PHASE PLANE AND TIME RESP9NSE 
PL9TS eN THE PRINTER. IN ADDITI0N ALL OF THE SYSTEM AND S9LUTI0N 
PARAMETERS are PRINTED S3 ONE KNOWS WHAT WENT INTO THE PLOTS. 
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21 IF(MeD(I/10).En»C)WRITE(l/20)L/T>(Z(N)/M=liN) 
IF(F • N E ♦ 0 • 0 ) WRUEU/20)L/Ti (Z(M)/M=1,N) 

0 FSRNAT ( I4#9F15*5 ) 

IF( 400- I >15,15/16 



15 WR I TE ( 6/ 10 ) T 

ENCODE (96/ 10, I TEXT ( 1 / 34 ) ) T 

CALL T E XTO ( I DEV/ I TEXT ( 1 / 34 ) / 24/ 34/ 1/ 1/ 2/ I ER } 

10 FORMAT ( ' SSLUTI9N TERMINATED AT TIME ='/F10.5/ 
*' WITH 400 SOLUTION POINTS') 
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